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Abstract 

The ground state of cerium metal in the intermediate pressure regime (5Gpa<P<12Gpa) is 
studied in detail by the newly developed LDA+Gutzwiller method, which can include the strong 
correlation effect among the 4/ electrons in cerium metal properly. Our numerical results show 
that the a" phase, which has the distorted body centered tetragonal structure, is more stable in 
the intermediate pressure regime and all the other phases include the a 1 phase(a-U structure), 
a phase (fee structure) and bet phases are either meta stable or unstable. Our results are quite 
consistent with the most recent experimental data. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.20.Bc 
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I. INTRODUCTION 



Due to the strong correlation effect, the 4/ electrons in Lanthanide metal usually par- 
ticipate very weakly into the chemical bonding, which makes these materials to be approx- 
imately s-band metal with close packed crystal structure. An very important exception of 
this qualitative understanding is the cerium metal, where the 4/ electrons do participate in 
chemical bonding in the a fee phase under ambient pressure. While the a fee phase is quite 
close to the instability, an isostructure phase transition happens by raising the temperature 
above 116K, after which the crystal structure remains unchanged while the volume expands 
by 16% and the 4/ e.ectrons become localisedfl. Father numerical studies by imp.ementing 
the first principle methods with dynamical mean field theory[2|, y[ show that the 7 phase 
may be stabilized by the entropy (4]. 

Another mysterious phenomena in cerium is the intermediate pressure phase. At zero 
temperature, the cerium metal form the face centered cubic (fee) structure for pressure below 
5 GPa and body-centered tetragonal (bet) structure for pressure above 12 GPa. But the 
experimental results of stable phases in the intermediate pressure region between 5 to 12 
GPa are still quite controversial {5- 10 1, as will discussed below in detail. 

Ellinger and Zachariasen applied X-ray diffraction studies on high pressure cerium with 
a diamond-anvil cellQ], and they reported that for pressure between 5 and 12 GPa, the 
orthorhombic a'-Ce phase with a-uranium type of structure, is the most stable phase. They 
also found another meta stable phase a" phase, which is monoclinic body centered with a 
deformed cubic face-centered structure. The conclusion that a'-Ce is the most stable phase 
between 5 and 12 GPa while a" phase is meta stable has been supported by some of the 
follow up experiments, i.e. jfj and [7|. While another group of experiments led to opposite 
conclusion, which indicated that the a" rather than the a' phase is the most stable phase in 
the intermediate pressure regime. Using diamond anvil cell and synchrotron radiation Olsen 
et al. studied the high-pressure phase diagram of cerium up to 46 GPa[8|. They reported 
that a" phase is the stable intermediate pressure phase of Ce, and no evidence for a' phase 
with a-uranium structure was found. This was the first experiment reporting the a" phase 
to be the stable intermediate pressure phase for cerium. After that several other groups also 
reported the similar experimental results suppor ting a" phase to be the most stable phase 
of cerium for pressure between 5 and 12 GPa 
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The first principle calculation is a powerful theoretical tool to predict the ground state 
phases of solid. During the past two decades, many efforts have been made to reveal the true 
ground state phase for cerium under intermediate pressure by first principle calculations. 
The early results of linear muffin-tin orbital (LMTO) method found the a' phase to be stable 
in the intermediate pressure regime between the low pressure fee phase and high pressure 



bet phase [U|, Il2j, which was consistent with the early experiments. After that, P. Soderlind 
and Eriksson et al. applied the generalized-gradient approximation(GGA) based on full- 
potential linear muffin-tin orbital(FPLMTO) method to the same problem|13j] and found 
that a" phase was the most stable phase only in a small pressure interval. P. Ravindran 
et al. have made systematic electronic structure and total-energy studies on Ce and didn't 
find any stable intermediate phases between the low pressure fee phase to high pressure bet 



phase 



141 ] . According to their results, both the a' and a" phases are meta stable with a" 



phase being lower in energy After that, LDA or GGA type calculations by other groups 



using plane wave method 



151 ] or exact muffin-tin orbitals(EMTO) method 16] also got similar 



results, that both a' and a" phases are not stable and the true ground states in energy of 
cerium are a phase (low pressure) and bet phase (high pressure). 

Although the 4/ electrons in the a phase are delocalized, the strong repulsive interac- 
tion among them still modifies its electronic structure significantly. Due to the insufficient 
treatment of the correlation effects, the bonding strength of the a cerium has been over 
estimated by the LDA type calculations, which leads to smaller volume and larger bulk 
modules comparing with the experimental data. In the present paper, we apply the newly 
developed LDA+Gutzwiller method, which can satisfactorily treat the strong correlation 
effects in the 4/ shell, to determine the ground states of cerium under pressure. We first 
apply the above method to study the ground state properties of a cerium under the ambient 
pressure. Our results show that both the volume and bulk modules are improved dramati- 
cally, which manifests the importance of the strongly correlation effect for the 4/ electrons 
in a cerium. Further we apply the same method to study the intermediate pressure phases 
of cerium, and the results show that in the intermediate pressure region the a" phase is the 
most stable phase and all other states are either meta stable or unstable, which is quite 
consistent with the recent experiments 

The rest of the paper is organized as following. A brief introduction of LDA+Gutzwiller 
method are given in Sec. [Ill In Sec. II III we discuss the main results of our LDA+Gutzwiller 
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calculations with the comparison to the recent experimental results and LDA(GGA) results. 
The summary and conclusions are given in the last section. 



II. LDA+GUTZWILLER METHOD 

Gutzwiller first introduced Gutzwiller variational approach to study the itinerant fer 



romagnetism in systems with partially filled d bands described by the Hubbard model 



17| 



Since then, the Gutzwiller variational approach has been widely applied to various of strongly 



correlated systems |18hl22l]. Recently, we developed a computational method to incorporate 



LDA with Gutzwil 
LDA+G hereafter) 



er variational approach named LDA+Gutzwiller method(simply called 
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"25||, by successfully applying to a number of typical correlated mate- 
rials, the reliability and feasibility of this method have been demonstrated. In the following 
we present the method briefly, please refer to our previous paper [25[ for more details. 

Similar with LDA+U or LDA+DMFT methods, in the LDA+G method, the LDA Hamil- 
tonian, which can be extracted from the first principle calculation, is implemented by a 
Hubbard-like local Coulomb interaction, which is not adequately treated within LDA. The 
effect of this local Coulomb interaction can thus be considered within the Gutzwiller varia- 
tional approach. The Hamiltonian can be usually expressed as 

H = H LDA + H int — H DC} (1) 

with 

Him = w f (2) 

where H^da is the LDA part Hamiltonian extracted from the standard LDA calculation, H int 
is the on-site interaction term, where a and j3 are combined spin-orbit indices of localized 
basis {<pi )Q } on site i, among which the local Hubbard interaction is implemented, a = 
1, . . . , 2N (N is orbital number, e.g. N = 7 for / electrons). H^c is the double counting 
term representing the average orbital independent interaction energy already included by 
LDA. Without the H int term, the ground state can be exactly given by the Kohn-Sham 
uncorrelated wave function(KSWF) l^o); which is a single Slater determinant made from 
the single particle wave functions. However, with the increment of the interaction strength, 
the KSWF is no longer a good approximation because it gives too much weighting factor for 
those energetically unfavorable configurations. In order to give a better description of the 



ground state, the weighting factor of those unfavorable configurations should be suppressed, 
which is the main idea of Gutzwiller wave functions(GWF) \^f G ). GWF is constructed by 
acting a many-particle projection operator on the uncorrelated KSWF, which reads 

\y G ) = p\yo) = Hmo), (3) 

i 

with 

Pi = ^ A ir m;r, (4) 
r 

m ir = \i,T) (z,T| , (5) 

where m^r is the projector to the specified configuration |T) on site i. In equation Q, the 
role of projection operator P is to adjust the weight of each atomic configuration through 
variational parameters Ajr(0 < A^r < 1)- The GWF falls back to KSWF if all A»r = 1- On 
the other hand, if A^r = 0, the configuration |T) on site i will be totally removed. In this 
way, both the itinerant behavior of uncorrelated wave functions and the localized behavior of 
atomic configurations can be described consistently, and the GWF can give a more accurate 
description of the correlated metallic systems than KSWF. 

The total energy of the above system can be expressed as the expectation value of the 
Hamiltonian equation ([1]) using GWF, which takes the form 

E To tai = (Vg\H\Vg) = (Vg\H L da\Vg) + (^ G \H in t\^ G ) - E DC , (6) 
In equation ([6]), the interaction energy is given as 

rm ir , (7) 

i,r 

where m^r is the weight of configuration F 

m iV = (^ G |m jr |^ G ) (8) 
According to equation ([3]), the LDA energy of equation ([6]) can be written as 

(y G \H LDA \V G ) = (y \PH LDA P\y ) = (*o|#lzmI*o>. (9) 
H-lda i s ca hed the effective Hamiltonian under Gutzwiller approximation. 
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The DFT calculations for realistic materials are always done in reciprocal space, so the 
formulas above should transform to the reciprocal space, We define the Bloch states of 
localized orbitals \ia) 

i 

Then H^ DA in k space can be written as 

Hlda = z a \ka) (ka\ + 1 — \ ka)(ka\) 

ka ka 

H LDA (Y,^\k'P)(k'P\ + 1-2 \ k 'P)( k 'P\) (11) 

k'p k'P 

+ - -2„)|A;a;)(A; / Q;|if LZ )A|^ / «)(^«|, 

kk'a 

where z a is renormalization factor for local orbital a which depends on those Gutzwiller 
variational parameters Ar, for those non-interacting orbitals, the corresponding z factor 
equals to 1. 

According to equation (0, EJ [H]), the total energy reads 

E To tai = (9f \(%2z a \ka)(ka\ + 1 - ^ \ka)(ka\) 

ka ka 

H LDA (Y^zp\k'/3)(k'/3\ + WI)l*o> (12) 

k'p vp 

+ X^ 1 _ Z D n « £a LDA + 2 Ermr _ 
a r 

where e^DA = J2k( ka \ H LDA\ka) and n a = ]Cfc(^o|&a) (ka\^ ). 

The total energy expressed in equation ffl2|) depends on both the uncorrelated " starting" 

wave function |\l/o) and those Gutzwiller variational parameters Ar, which can both be 

determined by minimizing the total energy. After we obtain the ground state wave function, 

we can calculate most of the ground state properties based on it, please refer to our paper 



251 ] for more details. 



III. RESULTS AND DISCUSSIONS 



Like the LDA+U and LDA+DMFT methods, in LDA+G method the Hubbard-like local 
Coulomb interaction U will be chosen as the only empirical parameter. For 4/ electrons in 
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cerium metal we set U=4eV, which is consistent with other LDA+U and DMFT calcula- 
tion. First we calculate the equilibrium volume and bulk modulus for a-Ce under ambient 
pressure to check the validity of the U value. Our results are shown in Tab. [I] together with 



the results from the all-electron FPLMTO calculations [26[], pseudopotential Plane- Waves 
calculations 15f, and experiment^, 8]. We can see from the table that the LDA calculation 
usually over estimate the bonding strength among cerium atoms, which makes the ground 
state volume obtained by LDA to be about 20% smaller and bulk modules to be 60% larger. 
After treating the correlation effect more carefully by Gutzwiller variational method, our 
results are in very good agreement with experiments. Thus we set U=4eV for all the cal- 
culations for cerium with different volume, and spin-orbital coupling effect is always fully 
included. 

Early experiments have reported two intermediate-pressure phases of cerium, a' phase 
and a" phase, together with the low pressure a phase and high pressure bet phase. The 
instability of these four phases under pressure is the main interest of the present work. 
Neglecting the tiny distortion in the a" phase, all these three phases including a, a" and 
the high pressure bet phase can all be treated within the same bet structure but with the 

nnn 

different ratio of the lattice constants c/a[8l. 191. Il3l|. which is illustrated in Fig. [TJ The c/a 
ratio is exactly \/2 for a phase with the fee structure and is found to be around 1.65 for the 
high pressure bet phase. The c/a ratio of a" phase is reported experimentally to be around 
1.5 < c/a < l.56\Q sj]. Therefore in the present paper, we first apply the LDA+G method 
to minimize the total energy of the cerium with bet structure respect to the c/a ratio as the 
function of volume, which mimics the competition among the a, a" and high pressure bet 
phases under pressure, and after that we will compare the total energy of these phases and 
the a' phase. Our main results have been plotted in Fig. [2] with the comparison to LDA 
and GGA. The results obtained by all the three methods agree quite well for volume smaller 
than 20A 3 indicating that the bet phase with c/a = 1.65(denoted by the green arrow in Fig. 
[2]) is the most stable phase, which is also quite consistent with the experiments 0, For 
volume between 20. 5A 3 to 22. 5A 3 , which corresponds to the intermediate pressure region, 
the results obtained by LDA, GGA and LDA+G are quite different. First of all, all the 
three methods predict that all the three phases appear as the locally stable phases in the 
intermediate pressure region. While LDA results indicate that the a" phase is only the meta 
stable phase in almost the whole intermediate pressure region except for volume being 2lA 3 . 
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And GGA can give a very small region, within which the a" phase is the most stable phase. 
Only LDA+G predicts that the a" phase is the most stable one of all the three phases in the 
whole intermediate pressure region. For volume larger than 23A 3 , the three methods again 
reach the same conclusion that the a phase is more stable and all the other meta stable 
states including a" and high pressure bctphases disappear completely, which is again quite 



consistent with the existing experiments 



io|. 



The optimized c/a ratio as the function of volume obtained from our LDA, GGA and 
LDA+G calculations are plotted in Fig. [3] together with the experimental results (§]. We 
can find clearly that the LDA+G calculation obtains the stable region of a"-Ce to be from 
20.5A 3 to 22.5A 3 , which is in good agreement with the experimental datajsj. While LDA 
and GGA only get a much narrow region for a"-Ce stable as shown in Fig. El 

The next issue to be addressed is the relative energy of the a' phase compared with the 
other three phases discussed in the previous paragraph. The a' phase has a orthorhombic 
a-U structure which can be viewed as distorted fee with some of the face-centered atoms 
being shifted from their original positions, as described by the parameter 2y 11]. The 2y 
value obtained experimentally by McMahon and NelmesjyJ is 0.2028A. If 2y=0.5 A and 
a=b=c, the standard fee structure can be restored. Therefore we can calculate the energy 
of a-Ce(fcc) and a'-Ce within the same a-U structure. We use alb = 0.5115, c/b = 0.8756 

n 

as obtained from the experiment |9|] and optimize the 2y value for any given volume. We 
find that 2y=0.21 gives the minimum-energy for volume being 22. 5A 3 , which is in good 
agreement with experiment joj] . We calculated the energy of a-Ce(fcc) and a'-Ce within the 
a-U structure frame for a volume region: 18 A 3 < volume < 28 A 3 . The energy difference 
between th.em(E a / — E a ) is plotted in Fig. @] together with the energy difference of E a n — E a 
and Ebct — E a obtained previously. Our results confirm that a' phase is always higher in 
energy. We thus rule out the a' structure as an stable intermediate pressure phase of cerium. 
This conclusion is in good agreement with the FPLMTO method calculations 14] . From Fig. 
H] we can also see that a" structure is the most stable phase among the above mentioned 
four possible phases within the volume region: 20. 5A 3 < volume < 22. 5A 3 . Therefore based 
on the LDA+G calculation, we conclude that a" phase is the most stable phase for cerium 
in the intermediate pressure region. This conclusion is quite consistent with the most recent 
experiments 0, 1(3 ] . 

In Fig. [5J we plot the renormalization factor of the 4/ bands in the a, a" and high 
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pressure bet phases as the function of volume. We can find that the renormalization factor 
of the 4/ bands decreases monotonically with the increment of the volume for all the three 
phases, which can be easily explained by the fact that increasing volume reduces the hopping 
integral between the neighboring / orbitals which enhances the correlation effect among 4/ 
electrons and thus reduces the corresponding renormalization factor. From the present 
LDA+G calculation, we find that the main consequence of the correlation effect in the total 
energy is to reduce the kinetic energy. Since the fee structure is close packed, compared with 
the a" phase, the a phase has relatively higher kinetic energy gain, which is over counted 
by LDA. Because of that, the reduction of kinetic energy gain captured by LDA+G is also 
more pronounced in the a phase, which raises the total energy of the a phase relative to the 
a" phase and makes it unstable in the intermediate pressure region. 



IV. CONCLUSIONS 



In summary, using the newly developed LDA+G technic, we have carried out systematical 
numerical study on the phase diagram of cerium metal under pressure. We found that 
the correlation effect among / electrons in cerium plays a crucial role to determine the 
most stable phase of cerium in the intermediate pressure region. The LDA calculation over 
estimates the chemical bonding contributed by the 4/ electrons, which leads to smaller 
equilibrium volume and larger bulk modulus comparing to the experimental data. With the 
increment of pressure, the overlap between 4/ orbitals becomes more and more pronounced, 
which reduces the correction to the total energy caused by the correlation effect. Therefore 
the correct description of the correlation effect which evolutes with the pressure becomes one 
of the key issues to obtain the correct phase diagram in the intermediate pressure region. 
Our numerical results obtained by the LDA+G method conclude that the a" phase is the 
most stable phase in the intermediate pressure region, which is quite consistent with the 
recent experiments. 
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FIG. 1: The relationship between the a phase (face-centered cubic, dashed lines with hatched 
circles— main cell), bet phase (body centered tetragonal, dash dot line subcell with unfilled circles) 
and the a" phase( C-face-centered monoclinic, solid line subcell with black circles). 
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FIG. 2: The total energy of body-centered tetragonal structure as a function of c/a axial ratio. 
LDA, GGA and LDA+G results are marked with black, red and blue line respectively, a, a" and 
bet phase are marked with purple, orange and green arrows. 
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TABLE I: Theoretical and experimental values of the equilibrium volume V and bulk modulus B 
for a-Ce from our LDA+G calculations, some LDA/GGA calculations 

r 

data 



20| and experimental 



LDA 

GGA 
Present (LDA+G) 
Experiment [5, 8] 



V(A 3 



22.74[26], 23.3[15] 
26.05(26], 26.3[15] 

28.91 

28-29 



B(GPa) 



60.5(26], 58.7(15] 
48.7(26], 43.0(15] 

36.2 

20-35 





1.70- 




1.65- 




1.60- 




1.55- 








1.50- 




1.45- 




1.40- 



— •— 


LDA+G 


— ▲ — 


LDA 


- T- 


GGA 


□ 


Expr.(1985)[8] 




17 18 19 20 21 22 23 24 25 26 27 28 



V(A 3 ) 



FIG. 3: The c/a axial ratio for the body-centered tetragonal structure as a function of volume for 
Ce. LDA+G results are marked with red solid line and filled circles, LDA results are marked with 
green solid line and filled triangle, GGA results are marked with blue dash line and filled triangle, 
experimental data (Ref. 8) are marked with black open squares. 
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FIG. 4: The relative total-energy curves with respect to the a phase(fcc). 
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FIG. 5: The average quasi-particle weight as a function of c/a axial ratio for angular momentum 
J=5/2 and J=7/2 /-electrons of Ce. 
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